Mesoscopic conductance fluctuations in dirty quantum dots with 



single channel leads 



Edward McCann and Igor V. Lerner 
School of Physics and Space Research, University of Birmingham, Birmingham B15 2TT, United 

Kingdom 

Abstract 

We consider a distribution of conductance fluctuations in quantum dots with 
single channel leads and continuous level spectra and we demonstrate that it 
has a distinctly non-Gaussian shape and strong dependence on time-reversal 
symmetry, in contrast to an almost Gaussian distribution of conductances 
in a disordered metallic sample connected to a reservoir by broad multi- 
channel leads. In the absence of time-reversal symmetry, our results obtained 
within the diagrammatic approach coincide with those derived within non- 
perturbative techniques. In addition, we show that the distribution has log- 
normal tails for weak disorder, similar to the case of broad leads, and that it 
becomes almost lognormal as the amount of disorder is increased towards the 
Anderson transition. 

PACS numbers: 72.15.-v, 73.20.Dx, 73.20.Fz 
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Recently it has been shown that the conductance of clean quantum dots with point-like 
external contacts (lead width w ~ hkp 1 , where hkp 1 is the Fermi wavelength) have non- 
Gaussian distribution functions JI],3- Weak transmission through the contacts means that 
the electrons typically spend more time in the system than that required to cross it so that 
the energy level broadening due to inelastic processes in the dot 7 <C E c where E c rs hD/L 2 
is the Thouless energy. This inequality corresponds to the zero mode regime which allows 
the use of non-perturbative techniques including random matrix theory || and the zero 
dimensional supersymmetric a model Q. In contrast, it is known that the distribution 
function of the conductance is mainly Gaussian || in a weakly disordered open sample 
connected to a reservoir by broad external contacts of width w > £, where I is the elastic 
mean free path. Inelastic scattering processes in the reservoir result in a level broadening 
7 ~ E c . 

In the present paper the aim is to determine the distribution function of a dirty quantum 
dot with two point contacts (which allow a single transport channel). We will describe 
the regime of continuous energy levels, 7 > A, where 7 is the broadening due to inelastic 
scattering in the dot and A is the mean level spacing. This overlaps with the supersymmetric 
(SUSY) calculations [[[]] in the ergodic regime, A < 7 < E c . While the SUSY approach is 
valid also in the quantum regime, 7 < A, a perturbative diagrammatic approach applied 
below can be used also for 7 > E c . In this case all diffusion modes contribute to the 
conductance rather than a single homogeneous "zero" mode which is the only mode taken 
into account within the nonperturbative SUSY calculations. The conductance distribution 
function has a non-Gaussian shape also for such a strong level broadening so that the non- 
Gaussian shape is due to geometric factors - namely, the point-like structure of contacts, 
rather than due to the dominance of the zero mode. In addition, we will use a standard 
renormalisation group technique to consider the role of increasing disorder in the dot. 

Traditionally the conductance of a system with broad, spatially homogeneous contacts is 
considered by means of the Kubo formula || . For a lead geometry which involves spatially 
inhomogeneous currents, however, the conductance is often more conveniently expressed via 
scattering probabilities using the Landauer-Biittiker formula 0. In this paper, we start by 
writing the conductance in terms of Green's functions with the help of the Landauer-Biittiker 
formula. Then we will determine the conductance distribution in the case of a continuous 
energy levels spectrum, 7 > A, finding the moments of conductance by diagrammatic per- 
turbation expansion in the parameter (7/A) -1 . Finally we will use an effective functional 
of the non-linear a model as a framework for the renormalisation group analysis necessary 
to describe dependence of the moments of conductance (and thus of the distribution) on 
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the disorder parameter. Following Ref. ||, we show that the nth order moments are pro- 
portional for large n to exp(un 2 ) (u is a certain parameter to be specified later) which is 
characteristic of a distribution function having lognormal tails. As the amount of disorder is 
increased towards the Anderson transition, the conductance distribution in the dot becomes 
almost entirely lognormal. Again, this is different from the conductance distribution in the 
ensemble of samples with broad leads which is also characterised by lognormal tails whose 
role is increasing with disorder, but remains mainly Gaussian within the whole range of 
validity of the renormalisation group analysis, even at the threshold of the transition. 

We consider weak coupling through the contacts from the disordered region to electron 
reservoirs and we label probabilities for tunnelling through the contacts as ati and cti which 
are assumed to be constant. Following [|] the level broadening due to inelastic scattering a 
(in units of the mean level spacing A) is chosen for convenience to be greater than ct\ and 
a.2- So we can write a = 27r 2 7/A where 7 is the total level broadening. We will present a 
perturbative calculation both in the many mode regime, 7 ^ E c , where the relevant small 
parameter is g^ 1 = (27r 2 E c /A)~ 1 (g is the average conductance of an open sample in units 
of e 2 /nh) and in the zero mode regime, A < 7 < E c , where the relevant small parameter is 
a -1 . The calculation is not valid in the region a < 1 which was the main area of interest for 
previous zero mode calculations Since we are modelling a disordered sample, and a ^ 1 
(the level spectrum is continuous), we do not include Coulomb blockade and electron-electron 
interaction effects. As a result, the calculations are applicable to a disordered sample whose 
electronic charging energy is negligible compared to the Thouless energy (i.e. with spatial 
dimensions larger than usual quantum dots). However the calculations are also applicable 
to quantum dots when the gate voltage is such that the addition of a single electron does 
not change the total energy, conduction occurs and disorder effects are relevant ||. Note 
that a similar non-Gaussian distribution of Coulomb blockade peak height fluctuations was 



found by Jalabert et al ||, and recent experiments |L0j appear to be consistent with this 
prediction. 

The framework for determination of the conductance is the Landauer-Buttiker formula 



We use it in the following form [] I 



G = |^£M + r*), (i) 



ab 



where the transmission coefficient is the probability of transmission from the channel 

labelled by a in the left (right) lead to the channel labelled by b in the right (left) lead. The 
conductance has been written explicitly in terms of transmission from the left and from the 
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right since this most symmetric form is required to consider the influence of broken time 
reversal symmetry on the fluctuations of the conductance. The transmission coefficients may 
be related to Green's functions by [11,12 

T% R) = -^G +H (ri,r 2 -,e)g~M(r 2 ,r i; e), (2) 

where Q + (G~) is a retarded (advanced) Green's function and ri, r2 are the positions of the 
point contacts. In the entire energy interval of interest, e — Ef ~ ^/t, the mean density of 
states Uq is a constant and the T^ R ^ are energy independent so we will subsequently drop 
the e label (r = £/vf is the elastic scattering time). 

Each point contact has a width w ~ hkp 1 which corresponds to a single channel only so 
that the point to point conductance, G, is obtained from the Landauer-Biittiker formula, Eq. 
(HI), with only one term in the summation. Ensemble averaged cumulants of the conductance, 
{{G n )), are given by 



e 2 ai« 2 



(( \Q + (ri, r 2 ) g~ (r 2 , r x ) + Q + (r 2 , n) Q~ (r 1; r 2 )l )). (3) 



h (hv ) 

We consider the point contacts to be separated by a distance greater than the mean free 
path so that |r x — r 2 | I. In Eq. (|3]) spin is explicitly included with an extra prefactor of 2 
and the term a,\oti represents the transmission probability through the contacts themselves. 
Ensemble averaging in Eq. @ is performed within the impurity diagram technique [13| 



and it is convenient to use a representation in which slow diffusion modes are explicitly 
separated from fast "ballistic" ones [FJ]. Fig. lb shows the dominant contribution to the 
mean conductance which contains one diffusion propagator (drawn as a wavy line which 
corresponds to a ladder series in the conventional technique flT3"|). At each end of the 



diffusion propagator there are two-sided 'petal' shapes which represent motion at ballistic 
scales since the average Green's functions {Q R ' A (r, r')) (drawn as edges of the petal) decay 
like exp(— |r — r'\/2£). At the diffusive scale, R ^ £, the petals reduce to the constant 
X2 = 27TUQT. The choice of diagrams is dictated by the inequality |ri — r 2 | 3> I. It means 
that a diagram with external points ri and r 2 connected solely by average Green's functions 
as in Fig. la is exponentially small. 

The diffusion propagator D(r, r'; to) can be represented at zero frequency as 

D ^ r < 0) = ? WTi exp|!q ' (r - r,)l s ^ aRl (4) 

where D = vprjd is the diffusion constant and, for a closed system, the summation is carried 
out over all q = nn/L, where n = (m, • • • , nd) are non-negative integers. In an open system 
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inelastic scattering occurs in the leads and the summation is cut off at low momenta q ~ L^ 1 
where L is the system size. For the system with point contacts the energy level broadening, 
7, due to inelastic processes inside the dot is inserted 'by hands' into Eq. (f|). As a result, 
in the many mode regime, 7 > E c , the summation in Eq. (^) may be approximated by an 
integration with cutoff at q ~ where the inelastic scattering length L in = (.D/7) 1 / 2 . 
This leads to 

2 / 2 W 2 , , (g^HUjR), d = 2 

C( R = 77-w5— n lffH MW« • ( 5a ) 

{Anyi^D \RL in J ^ (7r/2 ^ o) (£/jR)5 d = 3 

where i? = |r — r'|, g ~ {^ftY" 1 is the dimensionless conductance of an open cube of size £, 
d = 2 + e is the dimensionality of the dot, and K a is the modified Bessel function of the third 
kind of order a. In the zero mode regime, 7 -C E c , the summation in Eq. (U) is dominated 
by the q = term so that 

27T 

C(R) « — a < ,E C /A, (5b) 
a 

which means that the diffusion propagator is independent of R, spatial dimensionality, and 
the degree of disorder. 

Now the leading diagram for the mean conductance, shown in Fig. lb, can be evaluated. 
Substituting \2 — ^wuqt for the petals and ({R) for the diffusion propagator, one has 

(G) = ja ia2 C{R), (6) 

where R 3> £ is the separation of the point contacts. The mean conductance is proportional 
to C(-R) an d Eq. (j5b]) shows that in the zero mode regime, 7 < E c , the mean conductance 
depends on the level broadening 7, but not on the separation of the point contacts, the 
dimensionality, or the degree of disorder. On the contrary, in the many mode regime, 
7 > E c , the mean conductance depends on all these parameters via Eq. (^a|); since it is 
inversely proportional to go it actually increases as the amount of disorder increases. 

In order to calculate the variance of the conductance, we consider the following correlation 
function between the transmission coefficient from the single channel at ri to r2 and the 
transmission coefficient from the single channel at iV to r^', 



Ar(An,Ar 2 ) = 2 



e 2 aia 2 
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((£+ (r lf r 2 ) QT (r 2 , n) G+ (r/, r 2 ') QT (r 2 ', n0». (7) 



h (2tti/ ) 2 . 

We use the notation Ari = n — ri', Ar 2 = r 2 — r 2 ', where |Ari|, |Ar 2 | £. The main 
contribution to K is shown in Fig. 2a. It has two square "Hikami" boxes, Xii representing 
motion at ballistic scales which are connected by two diffusion propagators, and it gives 
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K(Ar 1 , Ar 2 ) =2 
X4(Ar) = 



h (2t(Uq) 



X 4 (Ar a ) X 4(Ar 2 ) [D (n - r 2 )] 5 



6;+ (k)£T (k) e* 



k.Ar 



(8) 
(9) 



To determine the behaviour of if on length scales |Ari|, |Ar 2 | <C £ we need to evaluate 
X 4 (Ar) accurately, rather than substituting it by const x5(Ar) which is sufficient for |Ar| ^> 
I. We find that 



J g + (k) Q- (k) e* Ar dk 



7T 



6/2 



TTZ/qT 



2k„Ar 



Hfl{k F Ar+iAr/2£) + Hfi{k F Ar-iAr/2i) 



(10) 



where i?h 2 are Hankel functions. In three dimensions this gives 



l e/2 



K(An,Ar 2 ) 



i2 r 



h 



aia 2 C,{R) 



sin 2 (fcpAri) sin 2 (fc^Ar 2 ) 



-(A ri +Ar 2 )/e 



11) 



(k F A ri ) 2 (k F Ar 2 ) 2 

0, to the correlation function for optical speckle 



which corresponds, in the limit |Ari| 
patterns found in flI5| . 

We find from Eq. (|iy) that X4(0) = {2t^vqt) 2 so that, in the single channel limit |Ari| = 
|Ar 2 | = 0, we get K(0,0) = (G) 2 /2. This corresponds to redrawing the diagram with the 
external points iV and r 2 ' exactly equal to ri and r 2 , respectively, as shown in Fig. 2b. 
Note that although in each of the boxes in Fig. 2a r m r' with accuracy up to £ (as the 
Green's functions represented by the edges of boxes exponentially decrease at scale £), such 
an accuracy would be insufficient for calculating the variance of the conductance of the 
point-contact dot as the area of order £ d ~ l would include j >1 channels. 

The variance is found by ensemble averaging Eq. ([3D for n = 2. Since there is symmetry 
arising from an overall exchange of spatial labels in any ensemble average, there are only 
two distinct contributions to the variance arising from the expansion of Eq. (|3]). The first 
is equal to the correlation function K(0, 0) from Eq. (0) contributed by 



term, [T^] 2 + 

the diagram in Fig. 2b containing two diffusion propagators. The second term, 2T^ b T^, is 
contributed by a similar diagram containing two Cooperon propagators. For broken time- 
reversal symmetry Cooperon diagrams are absent so that, overall, we get 



((G 2 )) 



(12) 



where (G) is given in Eq. (|B|). The factor (3 in Eq. (IT3) corresponds to Dyson's orthogonal, 
unitary, and symplectic ensembles: /3 — 1 in the presence of potential scattering only, (3 = 2 
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in the presence of a finite magnetic field that breaks time-reversal symmetry, and (3 = 4 in 
the presence of weak spin-orbit scattering. For (3 = 1 the result Eq. fll~2] ) is the well known 
large intensity fluctuations in speckle patterns |15|] . 

In order to determine the distribution function, we need a general expression for the 
main contribution to the nth cumulant. The leading diagrams are a generalisation of those 
for the variance, n = 2. For n = 4, for example, Fig. 3a shows a diagram contributing to 
a correlation function between transmission coefficients from four different input channels 
to four different output channels. It consists of two eight-sided Hikami boxes. In contrast, 
Fig. 3b shows a contribution to the fourth cumulant of the single channel conductance 
which has two 'daisy' vertices where each daisy consists of four petals. Similarly the leading 
diagrams for the nth cumulant are a generalisation of Fig. 3b with two daisy vertices, where 
each daisy consists of n petals, which are connected by n diffusion or Cooperon propagators. 
Each diagram gives a contribution of (G) n and a factor of (n — 1)! arises because of different 
ways of ordering the n propagators. For /3 = 1 all the diagrams from the expansion of Eq. 
(0) are present. However for (3 = 2 only the diagrams without Cooperons remain i.e. those 
that arise from [T^ b ] n or [T^] n . There are in fact only two such terms in the expansion for 
all values of n, whereas the total number of terms is 2 n . So the relative number of (3 = 2 
terms is 2/2™ = (3~^ n ~ 1 \ As a result, the main contribution to the nth cumulant is 

((G n )) = ^T^r- (13) 
This expression corresponds to the following distribution function 



G 0- 



i 



f(G) ' 



f3G 
W) 



(14) 



which is drawn in Fig. 4. For (3 = 1 the distribution peaks at zero conductance, whereas 
for (3 = 2 it peaks at (G)/2. It has the same form as the distribution for level width 
fluctuations of quantum dots in the resonance regime which was found in || based upon the 
hypothesis that chaotic dynamics in the dot are described by random-matrix theory. The 
result obtained here is based upon entirely microscopic calculations. However, for (3 = 2, it 
disagrees with the result of microscopic calculations by Prigodin et al PI within the SUSY 
approach. Their result for /3 = 2 is the same as our (3 = 1 result. This discrepancy 
arises from a different original definition of the conductance. Had we defined cumulants 
as averages of [T^ b ] n only, we would have the same result as in Ref. [[!]], as one expects 
in the region where both the exact zero-mode integration within the SUSY approach and 
straightforward diagrammatics are equally applicable. However the conductance is defined 
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TT| , |T2H as the sum of T^ b and T^, Eq. (JJ). When time- reversal invariance is broken by a 



magnetic field (i.e. for the (5 = 2 symmetry class), the left and right transmission coefficients 
are no longer equal for a generic asymmetric dot. Thus cross- terms like [T L ] m [T i? ] n " m no 
longer contribute to the nth moment of the conductance, producing the result different from 
the (3 = 1 case. It means that breaking time- reversal invariance suppresses small amplitudes 
in the distribution ( |T4D and increases the mean amplitude. This has already been noted by 
Jalabert et al || and we refer to their paper for further discussion. 

The distribution (14 ) is very simple but profoundly different from the conductance dis- 
tribution of an open system (with broad multi-channel external contacts). In the latter case, 
the variance is universal (of order e 2 /h) [[16,17] , and higher moments are much smaller than 
the variance so that the distribution is almost Gaussian ||. The tails of this distribution 
decrease, however, much slower r than Gaussian tails. We will show that this is also the 
case for the single-channel conductance distribution considered here. It is known || that 
expressions for cumulants of the conductance of an open system found in the lowest order 
of perturbation theory are not valid for n 2 > (q 1 where ( is the standard weak-localisation 
parameter: ( = ((R = £), Eq. fl5a|) . The reason is that the number of additional diagrams 
containing closed diffusion loops which describe higher order (in £ ) contributions to the 
nth cumulant increase so fast that it is n 2 ( rather than ( which takes the place of the 
effective perturbation parameter. We have found that corrections in powers of (o also arise 
in the present case of the conductance fluctuations of a system with single channel contacts. 
For example one such correction which consists of diffusion propagators and contributes to 
the 2T^ b T^ b term of the variance is shown in Fig. 5. Three similar corrections containing 
Cooperon propagators also occur so that, for the vertex corrections in the first power of £o, 
we get {{G 2 )) = (4/ j3 2 )(G) 2 ( . Similarly for the nth cumulant extra impurity ladders can 
be placed in n(n — 1) different places so that the corrections give a series of terms in n 2 £o, 
not Co- At large enough n, this enhancement of corrections by a factor n 2 means that the 
"main" contributions no longer dominate. 

In order to find expressions for large n cumulants we need to sum all the corrections in 
powers of (o which is not practical within the diagram technique. Instead the summation is 
performed using the renormalisation group procedure which is carried out in the framework 
of an effective field theory, a non-linear a model ||18|| , where averaging over realisations of 
disorder and averaging over fast degrees of freedom are performed in the derivation of the 
model. The averaging produces expressions for the nth cumulant of the point contact con- 
ductance in terms of functional derivatives with respect to a source field h(r) (for notations 
see §), 
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«n<*» = f2e2 aia2 r2 



v h (27rz/ ) 2 8iVV 
where (Zfh]) is a generating functional, 



ntr 



<5 2 



i=l 



v 5h i (ri)5h i (r 2 ) 



(Z[h]> , (15) 



M = i^:gi ' mh] = F[Q] + FJQ;h]. (16) 



Here the functional -F[Q] is a modification of the standard cr model functional, 
F[Q] = / d d r 



^Tr(VQ)Vr-^Tr(AQ) 
8 4 



(17) 



which takes account of the non-zero level broadening 7 (see discussion after Eq. (|])). The 
source field functional is 

F,[Q;h] = £ Fi m) [Q;h] = -° £ T m / Tr(hQrd rf r, (18) 

m=l ^~ m=l 

with bare values of the charges T m given by 

T (0) _ ( 2 ^- 3 ) !! (19) 
m ~ ml ' 1 j 

The Hermitian matrix Q obeys the constraints Q 2 = J, TrQ = 0. It may be represented as 
Q = Q^t m where r M are quarternion units and v = {AB; ij;pp'} stands for a set of additional 
matrix elements. The replica indices {AB} run from 1 to AT with the replica condition 
N = being applied to the final results, the loop indices {ij} label different conductances 
in the product Eq. (|I5D , and the indices {pp'} distinguish retarded and advanced Green's 
functions. These indices are required to eliminate terms in the perturbative expansion of 
Eq. flTSP which do not correspond to those in the standard diagram technique. The matrix 
source field h(r) is chosen to be Hermitian with the following pp' structure: 







®r + 


f ° 


h AB \ 




Wba 






\-h% A 







h = ®r + ®r 3 . (20) 



High gradient vertices g are not included in the functional Eq. flTT]): although they are 
involved in the renormalisation of the charges T m in Eq. ([19]) this could produce only a 
change in preexponential factors irrelevant here. 

The lowest order perturbational contribution to the nth cumulant arises from the term 
(F^IQ; h]) 2 in the expansion of Eq. (|i~6|). The vertex [Q; h] contains n source fields h(r) 
and thus corresponds to a Hikami box with n external points such as those in Fig. 2a for the 
variance and Fig. 3a for the fourth cumulant. This contribution is proportional to (T^ )) 2 
and does not reproduce the exact numerical coefficient of the diagram technique result, Eq. 
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(f[3|), for single channel contacts which arise from daisy vertices with a single external point 
only (Fig. 2b for the variance and Fig. 3b for the fourth cumulant). The reason is that 
the derivation of the a model involves averaging over fast degrees of freedom so that it is 
insensitive to details on local length scales (of the order of hkp 1 ). Nevertheless the a model 
accurately describes the behaviour of diffusive degrees of freedom which are the relevant 
ones for what follows. 

The renormalisation group procedure allows effective summation of the higher order per- 
turbative corrections which are logarithmic in 2d. The net effect is to substitute renormalised 
values of the charges for bare ones in the expressions obtained by the perturbative expan- 
sion of Eq. flTBD above. Results in higher dimensionalities can be qualitatively obtained by 
d = 2 + e expansion. 

The source field functional above, Eq. (|18D , is similar to the source field functional 
describing fluctuations of the density of states which is renormalised in ||. As a result of 
the renormalisation, the charges obey the following increase law, 



T n oc r£V (n2 - n) , (21) 



where 



M = ln^ = ln(l-Co) _1 (22) 
o 

In the weak disorder limit u ~ Co 1, whereas in the vicinity of the Anderson transition, 

f e In L/i, L < L c (a) 

u={ . (23) 

\\n{l- g c /g Q )-\ L > L c (b) 

Here a is the physical (renormalised) conductivity at length scale L and do is the classical 
(bare) conductivity at length scale i. L c is the correlation length which diverges as L c oc 
(fi'o — 9c)~ l ^ e i n the vicinity of the Anderson transition point go = g c . 

Substituting the renormalized charge in place of the bare charge in the leading pertur- 
bative results we get 

((G n )) ~ (G)™e 2u(n2 - n) , n > u~\ (24) 

This is valid for cumulants with n > u^ 1 whereas the universal expression, Eq. (|T3"|), is valid 
for n ^$ -ur 1 ! 2 . The exponential increase law for high n cumulants, Eq. (|24"D, is similar to 
that of the local density of states |fr9|,|5[ and it leads to lognormal tails of the distribution 
function 
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f(SG) ~ ^ ex P 



1 , 2 ( 
hr ' 



, ... „ . SG>{G)/u, (25) 

where 5G = G — (G). For weak disorder it ~ (o ^ 1 so the main part of the distribution 
is due to the low n cumulants and it has an exponential shape Eq. flU). Some very large n 



cumulants follow Eq. (|24| ) and the exponential distribution will have lognormal tails which 
appear for fluctuations SG ^ (G)/u. 

As the amount of disorder increases then u increases in magnitude, more of the cumulants 
follow Eq. fl24]) , and the lognormal tails become larger. Due to the condition of validity of 
the high cumulant expression, n > w -1 , the whole distribution will become lognormal in 
the region u ~ 1. This crossover from the exponential to the lognormal distribution occurs 
before the Anderson transition i.e. still within the metallic regime, since u = Ina^/a then 
u ~ 1 can occur for cr > a ^> 1. This is similar to local density of states fluctuations 
19| where a crossover from nearly Gaussian to completely lognormal occurs in the metallic 



regime for u ~ 1. 

Note that the lognormal distribution for local fluctuations originally obtained by the 
renormalisation group treatment has been rederived directly within the SUSY approach 
|20| . It is possible that the high gradient expansion [see note after Eq. (20)] corresponds 



to probing the new inhomogeneous vacuum found in ||20|| . However the new approach is 
applicable only to the weak disorder limit, u ~ Co, and could not describe the distribution 
of the many channel conductance. 

In summary, using diagrammatic perturbation expansion in the parameter (7/A) -1 , 7 > 
A, we reproduced the exponential distribution of conductance fluctuations in quantum dots 
with two single channel leads H in the zero mode regime, A < 7 < E c , and we demonstrated 
strong dependence on time-reversal symmetry. We have shown that the distribution has the 
same shape in the many mode regime, 7 ^ E c , but, in contrast to the zero mode regime, the 
mean and the variance are dependent on the spatial dimension, the degree of disorder, and 
the separation of the leads. Using the renormalisation group procedure we have shown that 
the exponential distribution has lognormal tails in both of the above regimes. As disorder 
increases, the lognormal asymptotics become more important and eventually there will be 
a crossover to a completely lognormal distribution. 
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FIG. 1. Lowest order perturbational contributions to the mean conductance: (a) exponentially 
small diagram for |ri — T2 \ 3> £, (b) dominant contribution. 




(a) (b) 

FIG. 3. Lowest order perturbational contributions to the fourth order correlation function, (a) 

multi-channel correlations; (b) fourth cumulant of single channel conductance. 
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FIG. 4. Point Contact Distribution Function. 
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FIG. 5. Lowest order in (q correction to the variance. 
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